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Abstract 


We study the classical and quantum mechanics of a two-dimensional ver¬ 
sion of a Fermi accelerator. The model consists of a free particle that collides 
elastically with the walls of a circular disk with the radius varying periodically 
in time. A complete quantum mechanical solution of the problem is possible 
for a specific choice of the time-periodic oscillating radius. The quasi-energy 
spectral properties of the model are obtained from direct evaluation of finite¬ 
dimensional approximations to the time evolution operator. As the scaled h 
is changed from large to small the statistics of the Quasienergy Eigenvalues 
[QEE] change from Poisson to circular orthogonal ensemble [COE). Differ¬ 
ent statistical tests are used to characterize this transition. The transition 
of the Quasienergy Eigenfunctions (QEE) is also studied using the test 
with n degrees of freedom. The Porter-Thomas distribution is shown to apply 
in the COE regime, while the Poisson regime does not fit the test with 
z/ = 0. We find that the Poisson regime is associated with exponentially lo¬ 
calized QEE whereas the eigenfunctions are extended in the COE regime. 
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To make a direct comparison between the classical and quantum solutions 
we change the representation of the model to one in which the boundary is 
fixed and the Hamiltonian acquires a quadratic term with a time-periodic fre¬ 
quency. We then carry out a successful comparison between specific classical 
phase space surface-of-section solutions and their corresponding quasi-energy 
eigenfunctions in the Husimi representation. 

I. INTRODUCTION 

In recent years a good deal of attention has been directed at trying to understand the 
quantum manifestations of classical chaos (QMCC). Although a complete understanding 
is not yet in sight, signihcant progress has been made in obtaining partial answers to this 
paradigm. This progress has primarily been achieved from studies of lower-dimensional 
models: 2-D for energy conserving models and 1-D for driven systems. In particular, very 
few studies have been carried out in two-dimensional time-dependent problems, for even in 
the classical limit the theoretical analysis is nontrivial. It is the purpose of this paper to 
consider a time-periodic two-dimensional problem which, because of the particular nature 
of the model, has a reduced degree of complexity and thus allows us to analyze its solutions 
in detail. 

The model considered here is a two-dimensional version of the thoroughly studied one¬ 
dimensional Fermi accelerator [|I|. The 1-D quantum problem was introduced in Ref. 0. 
There it was shown that for a very specihc form of the wall-oscillation function the Floquet 
evolution operator can be explicitly written and analyzed, without the complications of 
numerical time-ordering (which effectively precludes a thorough investigation). The model 
has been further analyzed by others The Fermi acceleration disk (FAD) studied here 

consists of a free particle bouncing elastically inside a circular boundary with a radius that 
oscillates periodically in time. Classically, the model exhibits a transition to chaos as the 
amplitude of the wall oscillation is increased. In the quantum case, the periodicity in the 
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wall oscillations allows us to use Floquet’s theory in terms of the one-period time evolution 
operator. 

This model is of interest for several reasons: hrst, because the dynamics is not kicked, 
as in almost all other time-periodic models considered, e.g. the thoroughly investigated 
periodically kicked rigid rotator model (PKR) [^]. Second, the explicit form of the tight- 
binding-like model obtained for the dynamical equations decays algebraically instead of 
exponentially, as in the PKR. The FAD model then allows for a comparison of our results 
to those of the kicked models, for it is important to know if the corresponding results 
are generic. Third, the possibility exists here that the spectrum of the evolution operator 
contains a continuous component [^], implying nonrecurrent behavior of the wave functions. 
Evidence suggestive of such a continuous spectrum has been reported in |Q. 

One of the major advances in the understanding of the QMCC comes from the clear 
differences found between the eigenvalue spectral properties of models that classically exhibit 
chaos and those that do not 0] . The tools used to measure these differences are extensions of 
those developed in Random Matrix Theory (RMT) [^]. In this paper we present an analysis 
of the spectral properties of the time evolution operator using these tools. 

While the FAD model provides considerable evidence for the spectral (or statistical) char¬ 
acterization of the transition from integrability to chaos, for reasons explained later, it turns 
out to be inconvenient in the study of the classical-quantum correspondence in phase-space. 
To overcome this difficulty, we map the FAD model onto a model of a particle subjected to 
a bounded, rotationally symmetric, inverted harmonic oscillator potential in the presence of 
time periodic kicks ( the Fermi inverted parametric oscillator or FIFO model) [Q. In this 
new representation, we can construct a mapping in terms of classical canonical variables, 
which in turn permits us to make explicit connections between classical phase space solutions 
and their corresponding quantum counterparts in terms of the Husimi distributions. This 
change in representation also allows us to make connections with experiments in mesoscopic 
quantum dots |^. 

The outline of the paper is as follows: In section II we dehne both the classical and 
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quantum FAD and FIFO models. The quantum problem is shown to be almost completely 
integrable for the specihc functional form of the oscillating radius. In Section III we present 
the bulk of our quantum results for the statistics of the QEE in the FAD representation, 
both for eigenvalues and eigenfunctions. Specihcally, we discuss the eigenvalue statistical 
properties in terms of different measures that include the nearest-neighbor spacing probabil¬ 
ity distribution and its integrated form, the A 3 statistic and the two-point function S 2 . In 
htting the data from the Poisson to the COE regime we use the Brody function as well as the 
Izrailev tests to try and quantitatively parametrize the transition region. We further discuss 
the eigenfunctions properties in different regimes using the distribution of v freedoms as 
a convenient parameterization of the results 0 , as well as their localization properties in 
terms of the time-averaged transition probabilities. In Section IV we discuss the explicit 
classical-quantum correspondences in phase space, carried out in the FIFO representation. 
We consider a series of classical orbits in the surface of section and then identify the Husimi 
contour plots of the eigenfunctions of the Floquet evolution operator that correspond to 
the classical solutions. Finally, in section V we present both our conclusions and also some 
questions left for future studies. 


II. THE CLASSICAL AND QUANTUM MODELS 

As we mentioned in the introduction we will use two representations of the model stud¬ 
ied here. Both representations have their own strengths and thus their analysis leads to 
complementary results. In this section we explicitly dehne the FAD and FIFO models in 
the classical and quantum limits. The classical mechanics of the FAD is discussed briefly 
while the classical solutions of the FIFO are discussed in more detail since it will be in the 
FIFO representation that the quantum-classical correspondence will be carried out. 
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A. Classical Fermi accelerating disk (FAD) 


The FAD model consists of a free particle confined to move inside a two-dimensional disk 
whose radius oscillates periodically in time; -R(r) = Rq6{t) with 6{t + T) = 6{t) and T is 
the period. The circular wall is taken to represent an infinite potential barrier for the motion 
of the particle, that is, classically the particle undergoes perfectly elastic collisions, while 
quantum mechanically the particle’s wavefunction vanishes identically at the boundary. 

The classical Hamiltonian written in polar coordinates is 
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( 1 ) 


where r is the radius with 0 < r < R{t), 0 the azimuthal angle, and m the mass of the 
particle. The resulting equations of motion are. 


(Pr 

dP 


m?r^ ‘ 


0 < r(r) < Rq9{t), 


(2a) 


d(f) 


J 


0 < 0(r) < 27r, 


(2b) 


dr mr‘^ ’ 

where J is the magnitude of the conserved orbital angular momentum. In principle, we 
can study the dynamics of the problem for any form of the wall oscillation function 0{t). 
However, as we shall show in the quantum case, the problem dictates the appropriate 9{t) 
that allows an explicit evaluation of the time evolution operator. The specihc form of 9{t) 
is found to be 


1 


l + 2e 

mod(r, -) 


( 3 ) 


where e is the amplitude of the wall oscillation. We will henceforth measure e in units of 
Rq and time in units of the period T, while the mass will be taken as one henceforth. Note 
that 9{t) is already discontinuous in its hrst derivative, which is important in order to have 
chaotic solutions. 

As mentioned above, the FAD representation is not the most appropriate one to discuss 
the classical-quantum correspondence of the solutions, of central interest in this paper. 
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However, to show that the FAD model does indeed show a transition to chaos we have solved 
the classical equations of motion at a period for different parameter values. Examples of Pr 
vs r phase space plots are shown in Figs, l(a-c) for three values of J and e. Note that since 
(j) is related to r(r) by Eq. (pF|), if r(r) shows irregular behavior, so will 0(t). 


B. Classical finite inverted parametric oscillator (FIFO) 


To obtain the FIFO representation of the model studied here we begin by carrying out 


the Liouville transformation 16 


T = j^e^{t')dt', 

r(r) =e{t)p{t). 

resulting in the new equation of motion for the p coordinate, 

(Po 

p -^ = 0, 0 < p{t) < Ro, 


(4) 

(5) 


where 


dR m?p^ 


d (Idd 


( 6 ) 


= -rA-n-n 


l(W 
6 dt 


dt \9 dt ^ 

The equation of motion can be derived from the transformed Hamiltonian 


(7) 


TT 1 2/ N 2 

= :^ + -mu, (t) p + —^, 


with 0 < p < Rq. 


( 8 ) 


This hamiltonian describes the one-dimensional motion of a particle in an oscillator potential 
whose frequency changes with time, as given by equation Eq.(|^), with a centrifugal barrier 
and conhned to a stationary inhnite well of width Rq. If, in addition, we choose the particular 
form of the wall oscillation function 6{t) as given in equation (j^), the explicit time-dependent 
frequency of the oscillator is given by 


OO 1 

u;\t) = -e^ + 26 ^ {5(f - (n + -)To) - 6{t - nTo)}, 

n=—oo ^ 


( 9 ) 


where 
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( 10 ) 


Thus, the oscillator frequency consists of a constant, negative part, and two sets of periodic 
Dirac ^-functions, one at integer, and the other at half-integer periods. This is the same 
form for Lj{t) obtained by Chu and Jose in their study the one-dimensional Fermi model [0. 
In normalized units, the full classical Hamiltonian Hpjpo reads 

„2 1 j2 oo 1 

HpiPO = — - - + ^ + 2ep2 £ - (n + -)) - - n)} . (11) 

The free parameters in the problem are e (associated with the amplitude of wall motion) and 
J (the angular momentum of the particle). Note that the harmonic oscillator kick strength 
alternates in sign every half-period. 

When the particle is free, that is when it does not hit the wall nor does it get kicked, 
the equation of motion is 


2 

p — e p-\ —j. 

The solution to Eq.(^) is easily verihed to be 
p{t) = 


\ 


Po 


Po cosh{e(t - to)} + -f sinh {e{t - to)} 


pit) = 


J 

+ sinh^ {e(t - to)}, 

Po^ 


( 12 ) 


( 13 ) 

( 14 ) 


or, in terms of the initial energy 


E« = 2 P? 




the solutions are given by 


pit) = 


\ 


a cosh 


pl+Eo/e^ 

2e(t — to) J- cosh (—^— 


pit) = ± 


\ 


V 

2Eq J- 


Eq 


(15) 


(16) 


(17) 
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The sign of p is positive or negative depending on whether the particle’s last collision was 
with the centrifngal barrier or the wall. Here, to, po and Eq are determined by the initial 
conditions, and the constant a is dehned as 






a = — + 


( 18 ) 


For a given energy E the radial coordinate is constrained by the wall and the centrifngal 
barrier: 


Pmin ^ P ^ 1 ; 


pmin 


\Ja - Ea/e^. 


( 19 ) 


Note that these eqnations of motion are non linear, even in the wall’s absence. Formally, 
one can represent either of the set of eqnations (^)-(14) or (|T6|) -(17) as 




pit) 
\Pit) j 


/ 

V 




( 20 ) 


/ 


E{po,po,t- to) 

G{po,po,t- to) 

bnt we can investigate the stability of the mapping by linearizing abont hxed points only 
nnmerically. 

The effect of collisions with the wall (or potential barrier) is simply to reverse the direction 
of motion of the particle : 


^ Pitt) ^ 


0 ^ 

' pK) '' 

^pitt) J 


o 

1 

^Pitc) y 


( 21 ) 


where tc is the time of collision of the particle with the wall (or potential barrier), and the 
minns and pins signs indicate times jnst before and after the collision. 

The effect of the kicks at every half period are obtained by integrating the eqnations of 
motion over an inhnitesimal kick time. The resnlts are, for half-integer and integer period 
kicks. 


Pitt) ] 
2 


1 ot 

%(tl) ^ 
2 

Pitp ^ 



(p(*T)) 


( 22 ) 
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( 23 ) 




1 oi 

^ 



V+2e ly 



respectively. The non-integrability of the problem arises precisely because of these kicks. 
In contrast to the one-dimensional (or J = 0) case, because the particle here can never be 
exactly at the origin at hnite energies, the kicks are always relevant (in the sense that they 
always couple the particle’s position and momentum). 

While so far we have only looked at the (p, p) plane, we can also take a coordinate-space 
cut. From the definition of angular momentum, 

d(j) J 
dt p^{t)' 

we obtain that during a ‘free’ propagation, the angle advances according to 

m -4 = 5 |sm-‘ (; (^ - f(; (yij5 - §) ) } (““d 24 (25) 

where a was previously dehned. In this paper we concentrate on the (p, p) phase space plane. 

We are now in a position to construct a map in phase space which takes an initial set 
of coordinates to a hnal one. The algorithm actually used was to propagate the initial 
conditions via Eqs. (pTD-(17) until it either hit a barrier, in which case Eq.(^) was used to 
change the phase space coordinates, or until the time elapsed was a half-integer or integer 
multiple of the period, in which case Eqs.(|2^) or (^) are used to transform the coordinates. 
We emphasize that in between kicks there may be many collisions between the particle and 
the walls, which effectively means new initial conditions for propagation under Eq. ([ibD - 
(17). In general, the map is quite complicated, and very sensitive to initial conditions. By 
recording the values at each successive period, we obtain a surface-of-section of the trajectory 
of the particle in phase space. 

We shall discuss specihc results under different conditions in Section IV. 
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C. Quantum FAD model 


The Schrodinger equation of the quantum mechanical version of the FAD (QFAD) model 
introduced in Section in polar coordinates is 


■ ^dA _ hM 02 Id 1 ^2 
dr 2m 1 dr"^ ^ ^2 g^jp. 


A(r,0,r), 


(26) 


with the boundary condition A (r(r) = Rq9[t),(I),t) = 0, and normalization 


RodA) 2n 

J rdr J d(j) |A(r, 0 ,r)|^ = 1 . (27) 

0 0 

These equations completely specify the QFAD model considered in this paper. To find the 
general solution to the QFAD model we expand the wave function T (r, 0, r) in terms of the 
natural complete set of basis functions ipni, 


H-C50 C50 

A(r,0,r)= ^ ^ C„i(r) V'„z(r,0,r), (28) 

l = — 00 71=1 

where the 'ipni^s are the instantaneous orthonormalized eigenfunctions for a given R{t), i.e. 
the solutions to 


h 

- 7 ^ V’ni (r, 0, r) = Eni (r) iJni {r, 0 , r). (29) 

2m 

Here r is a parameter and satisfies the same boundary and normalization conditions as 
A. The instantaneous solutions are given by 


0n0A 0, r) = exp(i/0), (30) 

y/TrR{T)Ji+i{l3nl) R{t) 

where Ji is the Bessel function of order I and jSni is its n-th zero. The corresponding 
instantaneous energies are 


Eniir) = 




Pnl 


2m \Ro6{t) 


(31) 


We now substitute the expansion (|^) into the Schrodinger equation (^ 6 [) . After using 
orthogonality and known identities involving Bessel functions, one can evaluate the integral 
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Ki{m,n) 


d 


= f Jl{l3mix) -^JliPnlx) dx 


(dml(dnl 

^ ^2 Jl—liPml) Jl—l{Pnl') 


Pnl ~ Pml 

SO that the final result can be cast in the form 


{m ^ n), 


ihCni = J2 


m,k 


9 2 (3^1^ Pnk 

Enl^k^l^n^m ^ ^7j ^k,l\^ C^nk 

Pmk ~ Pnk 


— S HjiYa-f^i Crnk'} 
m,k 


(32) 


(33) 


where the overdot denotes differentiation with respect to r. Equations (|33|) are completely 
equivalent to the Schrodinger equation (^). The time development of the system is obtained 
from 


Cnlir) = Unm-,kl{'^,T') Cmfc('r'), 

m^k 


(34) 


with the time evolution operator defined by 


Unm]kiij ~) Texp ^ ^ J^ ds^ , 


(35) 


where T formally represents the time-ordering operation. We can, however, get rid of the 
time-ordering if the effective Hamiltonian H factors into an overall time-dependent factor 
multiplying a time-independent part. For this to happen, we see from equations ( 0 ) and 
(B3|) that we must have 


9 1 

— rs_/ - 

9 02 - 


(36) 


An appropriate form for 9{t), satisfying the periodicity and time reversal invariance condi¬ 
tions is the one given in (§). Using this form for 9{t) the matrix elements of the one-period 
time evolution operator are exactly given by 


Uni-mk{r + T,t) = exp 


h 


H, 


(■’■+'^ ds 


nl]mk 


02 (s) 


(37) 
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Note that the time evolution matrix is nominally four-dimensional but, because of the con¬ 
servation of angular momentum, it reduces to a two-dimensional matrix, for a given value 
oil. 


The required zeros of the Bessel functions are calculated using the Newton-Raphson 
technique after bracketing the root using the fact that separation between successive zeros 
asymptotically approaches tt |T^. This method is very efficient, since hnding the roots even 
for the high zeros and large I values can be done to machine accuracy. The quasi-energy 
spectrum is determined once the zeros of the Bessel functions are found and the t/-matrix is 
diagonalized numerically (with unitarity preserved to machine precision). Since U is unitary, 
its eigenvalues are unimodular, and we write it’s eigenvalue equation as 


f/(r + T,r)|A„) = e-'*^|A„) 


0 < < 27r. 


(38) 


Complete knowledge of U implies full knowledge of the evolution of the system at all integer 
multiples of the period because U{t + NT,t) = U^{t + T, r). This product of f/’s is 
the quantum counterpart of the classical one-period maps. The quasienergy eigenvalues 
En, being invariant under unitary transformations, are representation independent and the 
quasi-energy spectrum obtained this way is of high quality from a numerical point of view. 
For our calculations, we used a truncated basis of 400 states for calculating properties of the 
QEF, while we used a basis of 500 states to calculate the QEE. We checked that the QEE 
spectrum was stationary, i.e. the eigenvalues remain basically unchanged to within 1% in 
this case as we increase the size of the U matrix from 400x400 to 600x600. 


D. Quantum FIFO model 

Quantum mechanically, the required set of Liouville and unitary transformations are a 
little more complicated (because of the additional requirement of normalization), but are 
nevertheless doable. We hrst transform the Schrodinger equation (p6D to an effectively one¬ 
dimensional problem by writing the wavefunction as 
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A(r,0,r) = 


$(r, r). 




vrr 


(39) 


the Schrodinger equation (^) now reads 




(40) 


We see from the above equation that the centrifugal potential barrier, (P — ^)?r/2mr‘^, 
diverges at the origin for I ^ 0. We hereafter consider only the nontrivial case of nonzero 
angular momentum, so we require in addition that <h vanishes at the origin. The <h normal¬ 
ization and boundary conditions are now 

rRo6(r) 


rmo(T) 

/ |$(r, r)|^ dr = 1, 

Jo 


(41) 


$(r = 0,r) = 0= ^{r = Rq9{t),t). 


(42) 


Equations (^0}]4^) dehne an effective one-dimensional problem for r. 

We now proceed to apply the Liouville transformations, Eq.(^. the transformed wave- 
function 4/(p, f) is obtained from the unitary transformation 0) 

(43) 

where the overdot denotes now differentiation with respect to t. Following the dehnitions 
and applying the rules of partial differentiation, we have 


$(r, r) = 


9{t) 


=exp i 


d dp d d d Id 

dr dr dp dr dt 9{t) dp 

d^ _ I d^ 
dr‘i 9‘^{t) dp‘^ 

d dp d dt d p9 d Id 
dr dr dp dr dt 9^ {t) dp 9‘^ {t) dt ’ 


(44) 

(45) 

(46) 


Applying these results to the set of equations (^0|- ^) , after some algebra, we hnally obtain 
Schrodinger equation satished by T(p, f). 


m 


d'^ 


2m dp^ 


ihJ^^!{p,t) = {-+ ]^mu‘^{t)p^ + 




2mp^ 




(47) 
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along with the transformed normalization and boundary conditions 


Jo 


(48) 


^ 0,^) ^ 0 ^ ^{p = R^,t). 


(49) 


The corresponding quantum FIFO (QFIPO) Flaniiltonian is 



(60) 


where 


j = - 1/4 


I > 0 . 


(51) 


The hamiltonian Hqfipq is simply the quantized version of the corresponding classical 
Hamiltonian. What is signihcant is that the frequency law for uj^{t) is the same in both 
classical and quantum mechanics. As pointed out earlier in [Q, this is a special effect arising 
out of the classical action consisting only of linear and quadratic terms. 

The one-period time evolution operator is then given (in units where the mass, time and 
length are scaled in units of m, Tq and Rq respectively) by 



(52) 


where 



(53) 


The hrst and third terms in Eq.([5^) correspond to the kicks at full and half periods, while 


the other two terms correspond to the evolution of the state under Hq. We note here that 
while the classical problem has two free parameters e and J, in the quantum case we have the 


additional parameter h. By writing the Schrodinger equation (^71) in terms of dimensionless 
quantities, it is easy to verify that the parameter dependence of U is via e, e/h and J. 

For convenience we choose a set of Fourier sine functions as the basis on which to diag¬ 
onalize the evolution operator to span the Hilbert space. While unitarity requires a large 
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number of Fourier components as the scaled h gets smaller, this choice of basis has several 
advantages -it not only ensures automatic satisfaction of the boundary conditions, but also 
considerably simplihes the numerical calculations involved. In a sense, the Fourier basis is 
also the natural choice to handle the very short-wavelength oscillations encountered in the 
QFAD model. The details of obtaining the matrix elements of the evolution operator in this 
basis are presented in the Appendix. 

III. STATISTICS OF QBE SPECTRA IN THE QFAD MODEL 

As mentioned in the introduction, one of the clear QMCC emerges when one compares 
the spectral properties of specihc model systems as appropriate parameters are tuned to 
classically produce a transition from integrable to completely chaotic regimes. In this section 
we follow the general thinking developed in RMT to implement different tests to quantify 
the spectral properties of the QFAD model. These properties are obtained from a direct 
diagonalization of the one-period time evolution matrix. For the results presented here we 
vary the value of J and e/h, while we take the wall amplitude hxed at e = 1. This value 
corresponds to a non-perturbative value, where the kinetic energy, the centrifugal barrier 
and the energy imparted by the wall all have about equal magnitudes. Furthermore, we 
rescale all lengths by e and all energies by h. The semiclassical limit of interest here is then 
obtained when / | cx), h J, 0, with J = \ kept constant. 

Next we discuss the RMT tests and their application to the results obtained for the QEE 
of the QFAD model. 


A. Nearest neighbor QEE distributions 

A local measure often used in RMT is the distribution of nearest-neighbor energy level 
separations, P{s), where s = Sn+i — £n- In the extreme integrable and chaotic regimes it 
has been established that P{s) takes the Poisson or Wigner distribution forms. 
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Pp(s) = e-‘, 

Pw(3) — 


TT --S^ 

f se 4 ^ , 


(54) 

( 66 ) 


respectively. A convenient and often successfnl parameterization of the P{s) obtained in the 


transition between Pp to Pw is provided by the Brody interpolation formnla |11|: 


Pa{s) = 7 (q( + 1) s“ exp(- 7 s"+^), 


(56) 


where 7 = [b (^) 


a+l 


, and r(x) is the Gamma fnnction. This distribntion is normalized 
and, by constrnction, has mean spacing (s) = 1. We recover the Poisson case taking a = 0 
and Wigner for a = 1. A criticism to the Brody distribntion is, however, that there is no 
first principles jnstihcation for its validity. The fact remains that it does fit the specihc 
resnlts found when considering explicit model systems. 

Instead of hrst looking at the P{s) we found that the hts are better if we start by 
calculating the integrated distribution 


n«(s) = [ ds'Pa{s), 

JO 


(67) 


where we use all the eigenvalue spacings calculated. The results for nQ,(s) are shown in Fig. 

-I—r CX+1 

2. Next we use the Brody form to find the exponent a to fit the data to Bq = 1 — e~'^^ 
The a exponent is determined from a linear ht to In \ln (1 — + 

The hts shown in Fig. 2 are obtained by using these values of a. These results clearly show 
the general trend. As the nonintegrability parameter e/h is increased, the n(s) goes from 
being Poisson-like to Wigner-like. For small values of s (< 0.1), the Poisson limit is well 
htted by the Brody function while for the COE the ht is not as good, as one would expect 
due to the level repulsion. The corresponding results for P{s) are shown in Fig. 3, with 
the Brody hts using the a values obtained from nQ,(s). We found that the hts are better 
by analyzing the data this way rather than hrst htting P{s) since in the binning process we 
lose information. 

An alternative phenomenological interpolation formula with some jnstihcation was pro¬ 
posed by Izrailev II who used the analogy between the partition function of a two- 
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dimensional Coulomb gas of charged particles on a circle (the Dyson gas) and the joint 
probability distribution of the quasi-energy eigenvalues of the COE to construct an approx¬ 
imate form for P{s) given by 


P,{s) = A - f) . 


( 68 ) 


where A{j3) and B{j3) are constants hxed from normalization and the condition (s) = 1. 
This formula also reduces to the appropriate Poisson and Wigner limits for /3 = 0 and 1, 
respectively. As we discuss below, this interpolation formula is also good for the QFAD but 
the Brody distribution provides a slightly better £t to our results. The interesting aspect of 
Izrailev’s distribution is that in principle it provides a quantitative connection between the 
energy level repulsion parameter (3 and the degree of localization of the QEF. The connec¬ 
tion is established by considering the ensemble averaged “entropy localization length” Is for 
the Un components of QEE for a chosen basis and dehned as (Is) = A^exp (^{Sn) — , 

where the information entropy 


N 


Sn{Ui, .., Un) = - X! In Wn, 


— \'Ur. 


(59) 


n=l 


Note that Sn is essentially the logarithm of the number of sites signihcantly populated by 
the corresponding eigenstate, and Sn‘^^ is the entropy corresponding to the random (COE) 
limit, introduced so that for completely chaotic states the exact limiting value of (Is) = N. 
The conjecture is that the quantity {ls)/N, which varies between 0 for completely localized 
states to 1 for fully extended ones, is identically equal to the repulsion parameter (3 in 

Eq.(ii. 

To check the validity of this conjecture, we htted the level spacing distribution to P/ 3 {s) 
given in Eq.(|58D and determined the value (3hist which minimized the We also calcu¬ 
lated the ensemble-averaged entropy localization length as dehned above, and calculated 
the corresponding f3s, and compared the two. The results are presented in Figure 4. We 
see that the agreement is good where the classical chaos is strong, but gets worse as the 
classical motion is more regular. Otherwise, the transition from COE to Poisson statistics 
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is consistent with the general trends observed previously. We conclude therefore that while 
the Izrailev model is intuitively more appealing, the evidence to support it is not compelling 
enough (in our model at least) to prefer it over the Brody formula. 


B. A3 and T? QBE statistics 


We also calculated higher-order eigenvalue spectral correlations. The average number 
of levels in an interval of length L is < n{L) >= where the <> stands for 

spectral average, and n{a, L) is the number of levels in an interval of length L starting at a 
and ending at a -|- L. Also important are the various moments of the level distribution. The 
one considered here is the second moment of the average number of levels in a given stretch 
of length L of the spectrum, the statistic 


< S^(L) >= (^{n{a, L)— < n{a, L) >Y'^ . 


( 60 ) 


Another often calculated statistic is the the Dyson-Mehta which measures the stiffness 

of the spectrum. This is dehned by 

ra+L „ 

A 3 (L, a) = yminA^B / [iV(a:) — Ax — B]^ dx, (61) 

where N{x) is the unfolded number density. In our case there is no need to unfold the 
spectrum. A 3 is just the least mean square deviation of N{x) from the mean straight line 
behavior. This statistic is directly proportional to the < > by Jq{L^ — sL‘^x+ 

x^)T?[x)dx, and thus can be calculated for the circular ensembles as well. The specihc 
theoretical predictions for the averaged < A^{L) > = are = 

■^(in{L} — 0.007, and < > = ^ . These results are correct in the asymptotic 

limit valid for 15 < L. 

In Fig. 5 we present our results for < A 3 > and < >. In these hgures one clearly sees 

the transition from Poisson-like (dashes) to COE-Wke (solid line) behavior as e/h is varied. 
We note that the A 3 statistic does not saturate in the COE limit, even for the maximum 
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interval L that we looked at, as would be expected from semi classical arguments originally 
proposed by Berry. Furthermore, for the largest L considered the Poisson limit does not 
present the knee seen in other completely integrable systems. All in all the results shown in 
Fig. 5 are consistent with what we have come to expect for the transition region. 


C. eigenfunctions distribution 


To study the behavior of the eigenfunctions we begin by discussing the statistics of 


the overlap of the eigenfunctions with the natural basis vectors. Several authors [|I^ have 
conjectured that as the classical motion changes from chaotic to regular, this distribution 
can be represented by a y^-distribution in u degrees of freedom, with u decreasing from 1 
(the Porter-Thomas limit) to 0 (the Poisson limit): 


/n\u/2 

PM = exp(-..j,/2), 


(62) 


where y =|< X\nl >p, |A) label the QEF and \nl) label a set of N orthogonal basis vectors. 
(The j/’s have been rescaled so that {y) = 1.) We have tested this hypothesis for the overlap 
strengths. The results, plotted on a logarithmic scale in Figure 6, show the general trend 
of decreasing v as the corresponding classical system becomes more regular. However, we 
note that as we move from the COE to the Poisson limit the hts to the get worse. Note 
especially the appearance of a sharp second peak well away from 0 as z/ decreases. This 
discrepancy is related to the fact that the results are basis dependent when not in the COE 
limit. Equivalently, we do not expect to have a good £t to the except perhaps if we 
take the special basis obtained from a semiclassical calculation. Even if we manage to get 
good agreement with the x^ for a properly chosen basis the result will not be generic and 
therefore the statistical analysis of amplitudes would lose its universal meaning. Universality 
does apply, however, in the COE limit. 
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D. Transition between localized to extended QEF 


The phenomenon of dynamic localization gives another QMCC |]^|. In this section we 
discuss the localization properties of the quasi-energy eigenfunctions when projected onto the 
natural basis. A physically meaningful quantity to calculate is the time averaged transition 


probabilities P{n,m) where {\n),n = label the basis states. While for small e/h the 

probabilities are peaked around the diagonal [i.e. the QEF are localized in energy space), 
for larger values the probabilities are spread out, indicating delocalized states. To quantify 
this transition we define the pair-correlation function C'(A) by 


1 N-A 

= + A = 0,l...iV-l (63) 

where P{no,m) = iimj.^ 00 ^ fi[ |< uo | m >p is the transition probability from the initial 
unperturbed state Iuq/) to the unperturbed state \ml). We find good fits of our evaluations 
of this correlation by the form. 


C(A) ~ 


(64) 


where is a localization length. The results are displayed in Fig. 7. We note that close 
to the near-integrable regime ^ << N, indicating localization of the overlaps, while for 
larger e/h, ^ > N. These results are consistent with the statistical analyses and with the 
expectation of a smooth transition from localized to extended solutions as we change from 
the Poisson to the COE regimes. We note that this results depends directly on th size of 
the U matrix being considered. 


IV. CLASSICAL TO QUANTUM PHASE SPACE CORRESPONDENCE 

In this section we consider the question of how to relate the quantum to the classical 
results. The natural habitat of nonintegrability in hamiltonian systems is phase space. As 
we mentioned in the introduction the FAD representation does not lend itself to an easy con¬ 
nection between the classical and quantum results [^. It is in the FIFO representation that 
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we can carry out this program. Here we start by discussing a set of explicit classical phase 
space FIFO solutions that we will then directly connect to the QFIPO model results. The 
connection is done by using the Husimi representation of the quasi-energy eigenfunctions. 


A. Classical FIFO model results 

In Figs. 8(a-c) we present results of the classical Poincare sections for three different 
angular momentum values, / = 1,5,10, for the parameter values indicated in the captions. 
We observe from these hgures that there is a symmetry line in the {p,p) section of phase 
space. As pointed out earlier [^] for the one-dimensional case, the reason for the existence 
of this symmetry line is as follows. Consider a particle kicked at f = 0. The position is 
unchanged, but the momentum changes to = p^~'> + 2epo- If we denote p^^'> by po) then 
at time 0*^“^ the particle had momentum = —2epo + Po- Since the Hamiltonian is time 
reversal invariant, we see that propagating a particle forward in time from {po,Po) is the 
same as propagating it backward from (po, 2epo —Po)- Thus, all results are symmetric about 
the line p = ep. This symmetry is also present in the quantum problem, as will be discussed 
below. 

Although the map is complicated, there are a few cases which can be studied analytically. 
The hrst of these is when the energy of the particle arises solely from its angular motion. 
The particle then does not feel the effects of the walls. The points po = epo lying on the 
symmetry line are then hxed points of order one of the map. To prove this, we begin by 
noting that the energy corresponding to this initial condition is given by Eq. ([T^) is 


En — 




.2 ’ 


( 66 ) 


while the momentum just after the kick at half period is related to that just before by 


p^i'^ = p[^ — 2e pi. 

9 9 2 


( 66 ) 


The energy just after the kick is then related to that just before by 
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(67) 

( 68 ) 




(pi * -2eptf - eVi + 

2^2 


Pi 

2 


= e[ ^ + 2epi (epi -p[ ^). 


The trajectory of the particle during the second half-period is simply the time reversed 
trajectory during the first half, so that the energy has to be the same before and after the 
kick. Thus, Eqs. (^) and (pT]) imply 


p[ ^ = epi 
2 2 


P^^^ = 
2 


■‘Pi¬ 


rn 


From Eq.(B3), and using the definition of momentum from Eq. 


Pi ^ = 
2 




2Eq + e^p\ 

2 pi 


we can deduce that 


(70) 



2 


( 71 ) 


Comparing Eqs.(|55[) and dUD we see that pi = po, and because of time reversal symmetry 
then tells us that at t = 1 we must have pi = po- The kick at full period now changes the 
momentum according to 


p’t’ =p[ ’ + 2 £po. 


(72) 


Using the same arguments for a full period as those used at half period, we can see that 



epo =Po- 


(73) 


Thus, (po, epo) is a fixed point of order one of the map, and numerically we can verify that 
it is of parabolic type, with a maximum allowable value of po for a given e and J. That 
is, this line of fixed points, which starts at {pmin, (-Pmin) ends at a point determined by the 
condition pi < 1. From substituting the initial conditions into Eq.(|I^) at half-period, we 
find that this point is given by 


max 

Po 


5-U2 


1 + 



(74) 
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The reality of this maximum imposes yet another restriction on the allowable values of the 
parameters, viz, 


-sinh(^)exp(e/2) < 1. 


(75) 


We now proceed to calculate one of the simplest nontrivial period-2 orbits, shown in Fig. 9. 

The particle starts at {po,Po) (shown as point A in the figure), with energy Eq. The initial 

conditions are such that the particle bounces off the wall at p = 1, reverses, hits the barrier 

at p = pmin (point C) and reverses momentum again before reaching the point (pi ,pi just 

2 2 

before the kick at half period (point D). Now, from Eq.|T^ we can write the time elapsed in 
motion under the influence of the potential only as 


f — to = ^ j cosh 


/p2 + Eo/e^ 


cosh 


-1 f Po + 


(76) 


Using this result, and substituting the appropriate p’s, we can write down the times elapsed 
during the three sections of the trajectory, Iab, Itscl = tcB and tcD- Equating the sum of 
these three times to 1/2, we can determine the position at D in terms of po and Eq. The 
result is. 


pi = \ a cosh{e -|- b{Eo)} 


En 


(77) 


where 


6 (Eo) = 2 < cosh 


-1 f Pmin + Eo/ 


e„s.rTi±^)Ucosh-('i±^). (78) 


Now, the relation between the momenta just before and after the kick, and the argument 
about conservation of energy (based on time reversal invariance) lead to the relation (O- 
Inserting this into Eq.(^), we hnally get 

En 


2En e 


+ — = a cosh{e -h b{Eo)}. 


(79) 


Since Eq is expressed in terms of (po,Po) and a = a(Eo), the above relation provides, in 
principle, a relation between po and po- 
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The particle, after the kick at t = 1/2, follows the time reversed trajectory to the point 
G (which is the same as A) at t = 1“ as shown, at which point it will have position po and 
momentum —po- Then it gets kicked, which now changes its energy to Eq. This time, the 
momentum after the kick is given by 

=-po + 2epo. (80) 

Using conservation of energy under time reversal again, we hnd that 

Eo = Eo + 2epo{epo-po), (81) 

which also changes the other functions of energy: 

a = vU/TUlUT. ( 82 ) 


Pn 


= Ua - Eo/e"^ 


(83) 


It is important to remember that both d and pmin are functions of po and po via Eq.(^). 
After getting kicked at H, the particle hits the barrier at I, reverses and gets to point J before 
getting kicked again, and following the time reversed trajectory back to the initial point A, 
as shown. The logic of determining the second relation between the initial conditions is the 
same as that used during the first quarter of its trajectory. By setting the total time taken 
\tHi\{= tin) + tij = 1/2, we hrst obtain a relation similar to Eg. (|77D , and repeating the 
argument about energy conservation due to time reversal invariance, we hnd 


En = 


2pj 


(84) 


with 


Pi^ ]ja cosh |e + c(£)o)} - 


(85) 


and 


c(Eo) = 2 cosh' 


-1 f Pmin + Eo/e'^ 


cosh 


-1 f Po + 


( 86 ) 
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Putting all these equations together, the hnal result is 

-^ + ^ = a cosh (e + c(Eo)|. (87) 

In principle, a physical solution to Eqs. (|79D and (^) may be found analytically for the hxed 
point {po,Po), in practice we hnd it numerically. The hxed point shown in the Fig. 8 (b) (for 
/ = 10, h = 0.026 and e = 1), labelled F, is given by po = 0.657385..., po = 2.659181.... 
When substituted back, this hxed point satishes the equations given above to machine 
accuracy. 

This analysis to determine this hxed point can be generalized to a family of hxed points 
of period 2 ■ [m + n + 1), with n, m = 0 , 1 ,... by letting the particle bounce between the 
wall and the centrifugal barrier (2m + 1) times during the hrst half period (between B and 
C), and {2n + 1) times during the third (between H and I). The equations to be solved are 
now 


j2 

+ — = a cosh{e + bm{Eo)}, 


2E, 


( 88 ) 


where 


bm{Eo) = (2m + 1) b{Eo) — 2m cosh 


-1 f Po + -^o/e 


(89) 


and. 


where 


F, 


2En e 


+ — = a cosh{e + c„(Fo)}, 


Cn{Eo) = (2n + 1) c{Eo) — 2n cosh 


-1 ( Pmin + -^o/ 


(90) 


(91) 


and Eq, a and pmin are obtained in analogy to the period-2 case. In principle, higher order 
hxed points (and their associated families) can be calculated analytically in the above manner 
by writing larger numbers of such coupled nonlinear equations. In practice, however, that is 
quite a formidable task. To determine the stable manifolds we exploit the existence of the 
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symmetry line around hyperbolic fixed points, which are simply the unstable manifolds under 
time reversal. In practice, numerical difficulties preclude determining the stable manifolds 
in any other manner. 


B. Correspondence 


We can now make a direct comparison between the classical and quantum FIFO results 
by employing a phase space approach. To do this, we use the Husimi representation of the 
QEF H. The Husimi distribution, interpreted as a probability density, is a coarse-grained 
version of the Wigner function which goes smoothly to the semiclassical limit. In practice, 
the most often used technique of coarse-graining is to take the overlap of the QEF with 
coherent oscillator states. For the radial coordinate the coherent state is 




„(P) = (^)i exp [-^y - p„)^ + f)} . 


(92) 


which is a minimum-uncertainty Gaussian wavepacket centered at (po^Fo), with root mean- 
squared deviations given by Ap = Y^h/2cr, Ap = ^Jfia/2, and a is the ‘squeezing’ parameter. 
This parameter is adjusted when making comparisons to the classical phase-space plots. The 
Husimi distribution of a single QEF 'ipe{p)) is then defined by 

.1 2 


•^£(Po,Po) = 


'^%,Poip)^e{p) dp 


(93) 


The Husimi distribution is obtained by scanning through the values of (po,Po) in the region 
of interest in phase space, and the result is compared with the classical surface-of-section. 
We begin the comparison by noting the symmetry about the line p = ep in the Husimi 
contour plots in Figs. lO(a-d). As mentioned earlier, this feature carries over from the 
classical results for the same reasons as there, and it is in fact used to effectively halve the 
numerical effort. 

All calculations reported here were carried out for the wall oscillation amplitude e = 1. 
In this case, all terms in the Hamiltonian are comparable in magnitude, which means that 
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we are in a non-perturbative regime. A few calculations were done for different e’s, but no 
new qualitative features emerged. In choosing a value of h, we were guided by the following 
considerations. The value of h has to be small enough so that the system is well into the 
semiclassical regime, yet large enough so that the dimension of the truncated Hilbert space 
N (which grows with decreasing h) is large enough to preserve unitarity. Moreover, N has 
to be such that the largest eigenenergy of Hpipo has to be larger than the maximum energy 
of the classical particle in the region of interest in phase space. These restrictions eventually 
led us to choose h = 0.026 for N = 100. All the interesting features seen in this model 
are manifested in this regime. Finally, the classical conserved angular momentum was kept 
identical to the quantum value, — 1/4. 

In general, during our numerical calculations, we found that about 80% of the QEE were 
‘reliable’, in the sense that their imaginary parts were < 10“^. This is one of the limitations 
of our numerical determination of the QEE, since being interested in the semiclassical limit 
we needed to choose small values of h (indicated in the right hand corner of Figs. (10-12). 
This meant that a large number of Fourier components were needed to accurately represent 
the QEE. The answers for the QEE were quite reliable when calculating the Husimi 
distributions. However, the result were not good enough to do a RMT statistical analysis of 
the QFIPO model spectrum. This is why we resorted to the QFAD model, already discussed 
in the first part of this paper. Since we had precisely the opposite difficulty with the QFAD 
model we see that these two representations lead to results which are complimentary. 

The classical analysis was carried out for different values of the angular momentum J. 
First, we iterated a single (arbitrarily chosen) initial condition several thousand times, which 
typically leads to the chaotic background as shown in the figures. Imbedded in this back¬ 
ground are KAM tori centered around elliptic fixed points, defined by choosing appropriate 
initial conditions. In Figure 2(a), we show several such tori, and in particular, a fixed point 
of period 2 (marked as FI and F2) which was determined earlier analytically. Also shown in 
each of the figures is a hyperbolic fixed point of order 6, marked by its stable and unstable 


manifolds. The fixed points were determined by using a modified Powell method of 
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determining zeros of coupled nonlinear sets of equations. This method, like all multidimen¬ 
sional root-hnding techniques, requires a good initial guess to converge to the hxed point, 
but once given that, determines the root and the Monodromy matrix (the Jacobian or the 
determinant of the linearized version of the map equations ) reliably and accurately. The 
hxed point is elliptic, parabolic or hyperbolic if the discriminant obtained from the eigen¬ 
values (he., {Traced — 4 • {Determinant)) is negative, zero or positive, respectively. In 
all cases, it was verihed within numerical error that the map was area-preserving, he., the 
determinant was equal to one. The unstable manifold was obtained by iterating the map 
along the direction given by the eigenvector corresponding to the eigenvalue larger than one. 
The stable manifold is given by the time reversed version of the unstable one. 

Comparison of the Husimi distributions tF'ip{po,Po) with the classical phase space plots 
show some striking similarities. There are, for many QEF, many structures which unmis¬ 
takably correspond to elliptic, parabolic and hyperbolic periodic orbits, as seen in Figs. 
10(a-c),ll(a-c) and 12(a-c). In the cases of / = 1,5,10, for example, the Husimi representa¬ 
tion of one of the QEF given in Fig. 12(a), 12(b) and 12(c) all sit on top of the analytic 
period-two hxed point marked as F. Also, seen in the hgures there are Husimis which peak 
exactly on top of the unstable hyperbolic period-6 hxed point, referred to in the literature 
as “scars”. This correspondence is so robust, in fact, that often when a good guess to 
the classical hyperbolic hxed points are unavailable, the Husimis are used as a guide to 
the location of the hxed point (being unstable, hyperbolic hxed points cannot be located 
without a very good initial guess). These enhanced probability densities are conjectured to 
play as important a role in quantum mechanics as the hyperbolic orbits play in classical 
chaos. Finally, a rare but persistent occurrence in all the cases considered is that of a single 
Husimi distributions peaked simultaneously over both elliptic and hyperbolic hxed points, 
rehecting a purely quantum-mechanical tunneling between the KAM tori. Here we have 
only shown representative results of the correspondence between Husimi distributions and 
classical solutions. 



V. DISCUSSION AND OUTLOOK 


In this paper we have presented results of a thorough analysis of the quantum manifesta¬ 
tions of classical chaos in a Fermi accelerating disk (FAD). We started by presenting results 
of the statistical properties of the quasi-energy spectrum for both the eigenvalues and the 
eigenfunctions. Most of the established tests of the manifestations of chaos in the quantum 
limit were implemented successfully with a few exceptions. To wit, we found clear evidence 
for the transition from Poissonian to COE-Wke behavior in the eigenvalue statistics when 
we changed Ti from large to small. We have also connected the Poisson to COE transition 
in the eigenvalues to the corresponding localized to extended regime of the eigenfunctions. 
In the COE regime the Porter-Thomas distribution gives a good account of the statistical 
properties of the probability amplitudes, while in the Poisson limit the test fails to give 
agreement with our results. 

While the statistical tests have become the hallmarks of the QMCC it is in phase space 
where a direct application of the correspondence principle can be probed. We found that 
the standard description in the FAD representation is not adequate to exhibit the quantum- 
classical correspondence. By carrying out a classical Liouville transformation plus the cor¬ 
responding quantum unitary transformation, leading to the Fermi inverted parametric oscil¬ 
lator (FIFO), the correspondence was implemented explicitly for a selected set of classical 
trajectories in terms of their Husimi quesienergy distributions in phase space. The classical 
dynamics was studied in terms of a stroboscopic phase space plot. Various periodic orbits 
were found either analytically or numerically for different values of angular momentum. The 
quantum quasi-energy eigenfunctions were represented in phase space using the Husimi dis¬ 
tribution. A direct correspondence was found between the classical stroboscopic phase space 
plot and the Husimi distributions of the quantum quasi-energy eigenfunctions. Periodic or¬ 
bits present in the classical phase space were identified with corresponding high probability 
distribution in the quantum quasi-energy eigenfunctions. We found the “scars” in the quan¬ 
tum quasi-energy eigenfunctions that correspond to the classical unstable periodic orbits. 
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We also found quantum quasi-energy eigenfunctions that corresponds to tunneling across 
KAM curves between stable and unstable periodic orbits. While classically the KAM curves 
block the chaotic diffusion, the quantum tunneling allows it. 

Our study extends the analysis carried out in one dimensional driven systems to quasi- 
two dimensional ones. This study is then the hrst step in the direction of testing the results 
hereby obtained with experiments, in particular in mesoscopic quantum dots in the presence 
of time varying magnetic helds. We shall discuss this problem elsewhere 
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APPENDIX: 

In this Appendix, we set up the matrix elements of the evolution operator U. The Fourier 
sine basis we choose is 


{p, 0|/, m) = \l — sin(m vr p) cos{l (p) 


(Al) 


The basis is orthonormal : 


{l',m\l,n) = 6u^i6m,r 


(A2) 


Then, the matrix elements of Hq are 


{l',m\Ho\l,n) = 


(A3) 


where 
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Tjmn 

-^0 


- “T 


1 

{m\ {p^ + e^p^)\n) H--(^ItI^) 

^ P 


^m,n 


^m..n 7 ' \ 


- 1 ) 


Amn 

m-n ^ ^ ^ 


TT 


(m^ — n^Y 

71^ 1 

I 4m^7r^ 


m = n 


{l^ - 

- 2 ^-^ {(^ + ^)Si ['7r(m + n)] — {m — n)Si [7r(m — n)]} 


where Si(a;) is the sine integral, 


s,M= r!EM*. 

Jo 


10 u 


We diagonalized Hq numerically to obtain its eigenvalues and eigenvectors 


Holfi) = hUf,\fj,) 

\p) = En=i{n\p)\n) 
= X]n=l Cn,fi\n). 


(A4) 

(AS) 


(A6) 


(A7) 


(A8) 

(A9) 

(AlO) 


Here we dropped the angular momentum quantum number index /, because angular momen¬ 
tum conservation forbids transitions between states of different 1. Using the completeness 
and ortho normality properties of the basis set, we obtain the matrix elements of U as 

N 

{m\U\n) = Y. x (All) 


E. 


p,q,r ^n,fi 


,{m\ \r){q\ |p). (A12) 


The last two quantities are again evaluated numerically by using methods for calculating 
integrals with rapidly oscillatory integrands. 
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FIGURES 

FIG. 1. Poincare surfaces of section in the {r,pr) plane, (a) Almost regular behavior for small 
wall amplitude (e = 0.001, / = 1). (b) Intermediate case where regular and chaotic motions coexist 
(e = 1.0,/ = 10). (c) Fully chaotic case (e = 10, / = 100). 

FIG. 2. (a)-(d) Transition from Poisson to Wigner for IIq as the tuning parameter is changed. 

The solid line in (a-d) corresponds to a ncj(s) fit with the a value obtained from the best straight 
line fit to £n£n [1 — 11^(5)]”^. The Brody fit is essentially on top of the calculated data. The values 
of were (a) (1.0, 10), (b) (0.5, 20), (c) (0.1, 100) and (d) (0.01, 1000) for e = 1. The repulsion 
parameter a for (a-d) are 0.10,0.20,0.60 and 0.86, respectively. 

FIG. 3. (a)-(d) Nearest-neighbor spacing distribution P{s) for the same parameters as in Fig. 
2. The dashed line in (a-b) corresponds to the Poisson limit, and to the COE case in (c) and (d). 
The solid line in (a-d) correspond to a Pa fit with the a values obtained in the fits of fig. 2. 

FIG. 4. (a)-(d) Fits of the data to the Izrailev formula P/ 3 {s) (solid line). The dots are the 
fits as given by (5s (see text for details). The best-fit and predicted values for {(5hist,(5s) are, (a) 
(0.20,0.05), (b) (0.22,0.12), (c) (0.48,0.50) and (d) (1.14,1.00). 

FIG. 5. The < A 3 (L) > (upper plot) and the S^(L) statistics (lower plot). The parameter 
values are the same as in Fig. 2 with (a) -|-, (b) □ , (c) o, and (d) for x. The dashed line is 
the exact Poisson result while the continuous line is the corresponding COE RMT results. We 
observe a continuous transition from Poisson to COE. In the Poisson case the calculated values 
are slightly above the predicted behavior but it does not show the bending characteristic of picket 
fence spectra. In the COE regime no knee is found either. 
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FIG. 6. (a)-(d) Distribution of amplitudes of QEF with the natural basis states for the same 
parameter values as in Fig. 2. Close to the COE limit (d) the amplitudes are nearly gaussian or 
Porter-Thomas randomly distributed. Away from this limit the distributions are not well fitted by 
the distribution with a significant difference seen close to the Poisson limit. This discrepancy 
is explained in the text. The values of u from the fits are, (a) 0.10, (b) 0.15, (c) 0.52, and (d) 1. 

FIG. 7. (a)-(d) Logarithm of the pair-correlation function defined in Eq.(64) as a function of 
A for the same parameter values as in Fig.2. The transition from localized to extended behavior, 
exemplihed by the value of the localization length ^ as compared to the matrix sizes (N=400 in 
this case). 

FIG. 8. Classical Poincare section for / = 1 (a), 5 (b) and 10 (c). The period-2 fixed points 
are denoted as F. The hyperbolic fixed points in (b) are obtained as explained in the text. Here, 
and in the following, h = 0.026, e = 1.0 

FIG. 9. Schematic representation of the trajectory of the particular analytically calculated 
period-2 orbit, as explained in the text. The vertical scale is in arbitrary energy units, and the 
small separations correspond to the same energy. The horizontal scale is in arbitrary length units. 

FIG. 10. (a-c) Contour plots of the Husimi distribution of a QEE which corresponds to the 

three classical solutions corresponding to Fig. 8(a) {I = 1). In (a) we have to Husimi distribution 
that corresponds to the period-2 solution calculated analytically in the classical limit, (b) Husimi 
distribution of another QEF which corresponds to the period-6 hyperbolic orbit marked by its 
stable and unstable manifolds - a ‘scarred’ eigenfunction, (c) A mixed case, where a single QEF 
corresponds to both a stable, period-2 elliptic orbit, and the unstable, hyperbolic period-6 orbit. 

FIG. 11. Same as in the previous figure corresponding to the case ^ = 5 of Fig. 8(b). 
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